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We consider cosmology in the Einstein- jEther theory (the generally covariant theory of gravitation 
coupled to a dynamical timelike Lorentz-violating vector field) with a linear jE-Lagrangian. The 
3 + 1 spacetime splitting approach is used to derive covariant and gauge invariant perturbation 
equations which are valid for a general class of Lagrangians. Restricting attention to the parameter 
space of these theories which is consistent with local gravity experiments, we show that there are 
tracking behaviors for the M field, both in the background cosmology and at linear perturbation 
level. The primordial power-spectrum of scalar perturbations in this model is shown to be the same 
that predicted by standard general relativity. However, the power-spectrum of tensor perturbation 
is different from that in general relativity, but has a smaller amplitude and so cannot be detected 
at present. We also study the implications for late-time cosmology and find that the evolution of 
photon and neutrino anisotropic stresses can source the M field perturbation during the radiation 
and matter dominated epochs, and as a result the CMB and matter power spectra are modified. 
However these effects are degenerate with respect to other cosmological parameters, such as neutrino 
masses and the bias parameter in the observed galaxy spectrum. 

PACS numbers: 04.50.-|-h, 98.80.Jk, 04.80.Cc 



I. INTRODUCTION 

For more than two decades Milgrom's modified New- 
tonian dynamics (MOND) [l| has been able to explain 
galaxy rotation curves which are conventionally consid- 
ered as an evidence of cold dark matter (CDM) on galac- 
tic scales. MOND modifies Newton's second law of mo- 
tion to /i(|a|/ao)a = — V$jVj where a and <i>7v are the 
acceleration and Newtonian gravitational potential, re- 
spectively; ii{x) is an effectively free function tending to 
unity in the limit \a\ 3> oq, with oq being a new fun- 
damental constant, which must have a numerical value 
of oo ^ (200 km/s)^/(10 kpc) in order to match obser- 
vations on a galactic scale. This theory looks like New- 
ton's when accelerations are large but is significantly dif- 
ferent when accelerations are small. On galactic scales, 
|a| ^ oq, so the Newtonian dynamics is modified, but in 
a way that can fit spiral galaxy rotation curves. Subse- 
quently, Bekenstein 0] built a relativistic theory which 
has MOND as a non-relativistic, weak-field limit, thus 
making the study of cosmology possible. In addition to 
the conventional tensor gravitational field, Bekenstein's 
theory involves a vector and a scalar field, and is there- 
fore dubbed TeVeS. Interestingly, it has been argued that 
TeVeS could also explain the large-scale structure forma- 
tion of the Universe without recurring to CDM [sj , thanks 
to the presence of the vector field 4] . 

Recently, the authors of Rcf. [5] showed that TeVeS is 
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equivalent to a vector-tensor theory of gravitation where 
the vector field has a non- fixed norm. They also showed 
that the correct MONDian limit could be realized with a 
single vector field having non-canonical kinetic terms . 
These results indicate that a vector field in the gravity 
sector might be an interesting component of the Universe 
(indeed the model in ^] and its generalized version we 
shall presented below in Eq. ([1]) could be used to explain 
dark energy and dark matter in background cosmology) 
and merits more detailed investigations. 

The idea of a vector field coupled to gravity has a long 
worldline (see for example Q for a review and [1] for fur- 
ther references), but in this work we will focus on the 
model described in which is the most well-studied 
one, and investigate its cosmological implications. This 
particular theory is based on a dynamical vector field cou- 
pled to gravitation that picks up a preferred frame and 
preserves general covariance. This vector field is unit- 
norm, timelike, and violates local Lorentz invariance. It 
is called the iEther field (or simply iE-field) and we will 
refer to the associated Einstein- iEther theory as iE-theory 
defined by the iE-Lagrangian. The iE-Lagrangian con- 
sidered in [9|] is a special case of our general model in- 
troduced in Eq. ([T]) below for a unit-norm vector field 
that includes terms up to second order in derivatives, 
and it has been extensively studied in various contexts 

In Refs. [l^ [l^, the background cosmology and pri- 
mordial power spectra of perturbations from inflation of 
a slight different model were also considered. Here, we 
investigate these for the model presented in Ref . [9| , and 
also study the evolution of linear perturbation to the M- 
field during the radiation and matter-dominated epochs. 
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As we will show below, if we restrict the parameter space 
of the underlying theory so as to satisfy the local exper- 
imental gravity constraints, [13, [13], this perturbation 
becomes sourceless and decays during the epoch of infla- 
tion and late-matter domination. However, it is sourced 
by the evolutions of the photon and neutrino anisotropic 
stresses during the radiation era and early matter era, 
which have some imprints on the cosmological observ- 
ables. 

Our presentation is organized as follows. In Sec. [Ill we 
briefly introduce the general iE-theory and derive per- 
turbation equations for the background Friedmann-likc 
cosmology in the covariant and gauge invariant (CGI) 
formalism (see [l3l for a derivation of the perturbation 
equations in conformal Newtonian gauge). In Sec. IIIII 
we shall use these equations to discuss the perturbation 
dynamics for the cosmological models of Ref. [^. First, 
we summarize the existing constraints on the model in 
Sec. nil A[ then, in Sec. IIIIBl we present the evolution 
equations for the perturbation variables and then use 
them to show how the primordial spectra of scalar and 
tensor perturbations in this theory are unmodified and 
modified, respectively, on comparing them with the pre- 
dictions of general relativity (GR) . The late-time evolu- 
tion of the ^E-field perturbation and its effects on cos- 
mological observables are also studied there. Finally, our 
discussion and conclusions are presented in Sec. IIVI 

Throughout this work our convention is [Vq, Vf,]?/*^ = 
^ab d'^''' ^ Rab = Racb^ whcrc Rabcd, Rab are respectively 
the Riemann tensor and Ricci tensor; the metric signa- 
ture is (-I-, —,—,—) and the universe is assumed to be spa- 
tially flat, filled with photons, baryons, CDM, 3 species 
of neutrinos and a cosmological constant. 

II. FIELD EQUATIONS OF EINSTEIN-iETHER 
THEORY 

In this section we briefly introduce the general form of 
the JE -theory and derive the CGI perturbation equations 
which we will use to analyze the cosmological effects of 
the JE -field. The equations presented here are for general 
Lagrangians. Later, we shall focus on a specific class of 
such theory characterized by a linear Lagrangian (with 
f{K) = ~K, see below). 



gravitational action of the form 
IdttGn J 

in which Gn is the bare Newtonian gravitation constant, 
A is a Lagrange multiplier ensuring that JE -field has a 
unit norm, and Cje is the M -Lagrangian expressed by 

Cm = /(if'^'',rfVa^"V6^'^) (2) 

with 

cd = cig^'gcd + C26^A + c^S^Sl + c^JE'g.d (3) 

and /(■••) being an arbitrary analytic function of its ar- 
guments. As long as the norm of JEa is fixed {JE°'JEa = 1), 
the form of if" ''^.^ is the most general possibility we can 
have for our vector-field Lagrangian. Notice that Eq. ([3|) 
here differs from that given in [^ by the C4 term and 
we shall refer to our model and that of [3] as the GEA 
(Generalized Einstein iEther) model to distinguish from 
the one considered in ^] (for which we instead called EA 
(Einstein iEther)). The matter Lagrangian is taken to be 
the same as in standard ACDM model. 

We treat the JE -field JE"' and inverse metric g"^ as the 
dynamical degrees of freedom and vary the action with 
respect to them to obtain the field equations. The former 
gives the JE -field equation of motion (EOM): 

\/b{FJ^^) - C4FJE''\/bJE'\/aJEc = XJEa (4) 

where we have defined K = K''''^jW aJE^^W bJE'^ , F = 
F{K) = df{K)/dK and J"^ = K^^^S/ bJE'^ ■ The vari- 
ation with respect to the metric leads to a modified Ein- 
stein equation. One could retain the form of Einstein 
equations in standard GR by treating the vector field as 
a new contribution (denoted by T^) to the total energy- 
momentum tensor in the universe, in addition to that of 
the conventional fluid matter which is denoted by T^^^. 
Then, according to the definition 



A. The General Einstein-jE ther Theory IQ-kG^T^ = - 

The model we consider here is a slight generalization 
of that presented in [6]. It is characterized by a general in which £^ = Cm + ^{JE^JEa — 1), we have 



JE'J, 



(ab) 



^{aJ% - M{aJb) 



} + Igabf + [jEdVciFr'') - C4F(^''Vd^e)(^'Ve^") 



-F 



CiV,^"Vb^e - CiV=^,V,^fc - C4(^"Ve^a)(^''Vd^fc) 



(5) 



3 



We also note that by varying the action with respect 
to the Lagrangian multipher A we simply get the normal- 
ization relation of the JE -field, JE°'JEa = 1, as mentioned 
above. 



B. The Perturbation Equations in General 
Relativity 

The CGI perturbation equations in general JE -theories 
are derived in this section using the method of 3 + 1 de- 
composition [Ts'l (see fl§\ for applications of this method 
in modified-gravity models). First, we briefly review the 
main ingredients of 3 -f 1 decomposition and their appli- 
cation to standard general relativity [l8| for ease of later 
reference. 

The main idea of 3-1- 1 decomposition is to make space- 
time splits of physical quantities with respect to the 4- 
velocity w° of an observer. The projection tensor hab is 
defined as hab = dab — UaUb and can be used to obtain 
covariant tensors perpendicular to u. For example, the 
covariant spatial derivative V of a tensor field T^.'.'.g is 
defined as 



k 



(6) 



The energy-momentum tensor and covariant derivative 
of the 4-velocity are decomposed respectively as 

Tab = T^ab + "iqiaUb) + pUaUb - phab, (7) 
VaUb = C^ab + ■^ab + ^Ohab + UaAb. (8) 

In the above, nab is the projected symmetric trace-free 
(PSTF) anisotropic stress, qa the heat flux vector, p the 
isotropic pressure, aab the PSTF shear tensor, ujab — 
V[(jMf,] the vorticity, 9 = V^Uc = 3d/a (a is the mean ex- 
pansion scale factor) the expansion scalar, and Ab = iib 
the acceleration; the overdot denotes time derivative ex- 
pressed as (/) = u°''SI a4>, brackets mean antisymmetrisa- 
tion, and parentheses symmetrization. The 4-velocity 
normalization is chosen to be u"Uo = 1- The quantities 
T^ab, Qa, P,P are referred to as dynamical quantities and 
<^ab,'^abT 9, Aa as kinematical quantities. Note that the 
dynamical quantities can be obtained from the energy- 
momentum tensor Tab through the relations 



p 


- Tabu''u\ 




p 






qa 








= KhfT^d+phab. 


(9) 



Decomposing the Riemann tensor and making use the 
Einstein equations, we obtain, after linearization, five 



constraint equations 



= ^^{e''\yvOab): 



nqa 

Bab 



h V (Jab + V nJab, 



b — 



V'Bab 



6, 



— K 

2 



2 



^b)ec^^'^ 



Vc(7d + {p + p)^cd 



cd^,b 



(10) 

(11) 

(12) 
(13) 
(14) 



and five propagation equations, 

0+lo^ -^^Aa+'^{p + 3p) = 0;(15) 
2 - 1 

<^ab + i:0(7ab - ^ {aAb) + Eab + -^KT^ab = 0;(16) 



w+-euj-V^aAb] = 0;(17) 







1 


T^ab H 


- ^d-Kab 









{p+p)<Jab + ^{aqb) 
tab + OSab ~ ^^Bdla£, 



^d{a^b)ec 



Bab + OBab + ^"£d(ae 



b)ec 



+ Trd{a£b)ec 



0;(18) 
0.(19) 



^u'^u'^Eac^yVefbd- The angle bracket means tak- 



Here, Sabcd is the covariant permutation tensor. Sab and 
Bab are respectively the electric and magnetic parts of 
the Weyl tensor Wabcd, defined by Sab = u'^u'^^Vacbd and 

Bab 

ing the trace-free part of a quantity. 

Besides the above equations, it is useful to express the 
projected Ricci scalar R into the hypersurfaces orthogo- 
nal to u° as 

R = 2kp- -9^. (20) 
3 

The spatial derivative of the projected Ricci scalar, r]a = 
iaVai?, is then given as 

Va = KVap-y^Va^, (21) 

and its propagation equation by 

29 2 
??a + y7?a = --9aVay ■A-anVa'^-q- (22) 

Finally, there are the conservation equations for the 
energy-momentum tensor: 

p+(p + p)(? + V>, = 0, (23) 

qa + ^9qa + {p+p)Aa-VaP + y''Trab = 0. (24) 

As we are considering a spatially-fiat universe, the spa- 
tial curvature must vanish on large scales and so i? = 0. 
Thus, from Eq. we obtain 



1 



-6*2 = Kp. 



(25) 
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This is the Friedmann equation in standard general rel- 
ativity, and the other background equations (the Ray- 
chaudhuri equation and the energy-conservation equa- 
tion) can be obtained by taking the zero-order parts of 
Eqs. (HailSl), yielding: 



3 2 

p + {p + p)e 



(26) 
(27) 



In what follows, we will only consider scalar pertur- 
bation modes, for which the vorticity zuab and magnetic 
part of Weyl tensor Bab are at most of second order [l^ , 
and so will be neglected in our first-order analysis. 



C. The Perturbation Quantities in JE -Theory 

In the Einstein-jEther theories where we consider the 
M -field as a new species of matter, the gravitational field 
equations Eqs. ((T0]-[27]) listed above preserve their forms, 
but the dynamical quantities p,p, qa, i^ab appearing there 
should be replaced by the effective total quantities of the 
same type. For simplicity, we shall always use variables 



without superscripts to denote these effective total quan- 
tities, while for those of a specified matter species we 
shall add corresponding superscripts {e.g., denotes 
the energy density of the M -field • • • ) . 

The vector JE -field, JEa, requires further discussion. 
As we mentioned above, it has the normalization rela- 
tion M°'JEa = u°'Ua = 1. In the background Friedmann- 
Robertson- Walker (FRW) universe the requirements of 
homogeneity and isotropy require that lEa is just equal 
to Mq, which is unambiguously chosen as the 4- velocity of 
the fundamental observers. But in a perturbed, almost- 
FRW, universe this is no longer true and we can write 
= Wa + aSa, — u"" + ic"" whcrc is another (first- 
order) vector field that vanishes in a FRW Universe: we 
call it the perturbation of the iE -field. Furthermore, 
the relation M^lRa — u'^Ua = 1 implies that u"aGa — 0, 
i.e., aSa is a spatial vector field which is perpendicular to 
Ua, up to first order in perturbation. This fact is used 
extensively in deriving the perturbation equations {e.g., 
V"a3a = V°a3a etc.). 

With these preliminaries at hand, and after some 
lengthy manipulations, the JE -field EOM Eq. ^ can 
be written as (up to first order) 



J 



Ci4 



F&a +[f + sea 



Ci4 



FAa + FAa + -FOAa 



^{a-cu){Fe + Fe)~'^cuFe'' 



^aVa{Fe) + ^aFVaV^Sib + CigFV'' (aab + V(a 



0, (28) 



r 



and from the definitions Eq. ([9]) the JE -field energy den- stress can be identified from Eq. fTl (again up to first 
sity, isotropic pressure, heat-flux vector and anisotropic order) as 



np^ = ^/-^Fa(6l^ + 26lV°a3a)-Hci4FV° (yla + aca + ^6'ffia 



up 



JE 
'ab 



1 



1 



1 



2/ + 3 ^(^ + ^"^a) + 3 F 9 + e-' + {V-aSa)- + 2eV''s,a 



= -Cl4 



FAa + FAa + ^FOAa 



Ci4 



-ci3{F + F9) 



(Jab + V^affib) 



Fa2a + iF + Fe]sCa 



^a\/a{F0) + iaFVaV''a3fc + C13FV'' ((Tab + ^ (aS^b) 



(Jab + {"^{a^b))' 
I 



^{a-cu){Fe + FB)-'^cuFe'' 



(29) 
(30) 
860(31) 

(32) 
(33) 



where, in Eq. (|32p. we have used Eq. (|28|) . Here, we have 
defined the new parameters a = ci+2tC2+c^, C13 = ci+c^ 
and Ci4 = Ci -I- C4; K is given by k = SttGn. Including 
these JEj -contributions to Eqs. ([10] -Ell) o'^s obtains the 



modified gravitational field equations for the general JE 
-theory. It is also easy to check that the above results sat- 
isfy (separately) the conservation of the JE -field's energy- 
momentum tensor, Eqs. ([23l IM]) . Note that Eqs. 
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I33p are the general expressions of energy density, pres- 
sure, heat flux and anisotropic stress in the 3-1-1 decom- 
position which include both zeroth order (background) 
and first order terms; to calculate the actually density 
contrast etc. (see Eqs. ((5l1 - [54)) below) one needs to take 
the covariant spatial derivatives of these equations • 



and in the latter environment the observed gravitational 
constant Goo is also a rescaling of Gn, but with a different 
factor fH 



Go 



1 



1 



-G]\ 



(36) 



III. A SPECIFIC MODEL: THE LINEAR 
LAGRANGIAN 

In above we presented the field equations for general 
M -theories, but in what follows we shall only analyze 
the cosmology of a specific edition of the theory which is 
defined by choice of a linear Lagrangian: 



fiK) 



-K 



(34) 



(note that the minus sign is because of our sign con- 
vention). This model is by far the most well known, in 
the sense that it has been investigated in the contexts 
of static weak-field limit [l3| (see also for the weak- 
field limit of the iE- model considered in [6[), background 
cosmology [l^ , the radiation and propagation of the iE - 
gravitational waves , [l^ , compact stars [ll] , and black 
holes [IJ] . Some of these studies have imposed stringent 
constraints on the viable parameter-space of c^s. In view 
of these restrictions we will confine our study to this con- 
strained subset of possible theories. Note that the per- 
turbation dynamics of the M -model has also been ana- 
lyzed in [l^ in the absence of the C4 term. Here, we shall 
include this term and perform a similar analysis but in 
slightly different manner and in more detail; we will also 
discuss on some detailed features which lead to our con- 
clusions being different. The late-time evolution of the IS, 
-perturbation is also investigated. In particular, we shall 
find that within the locally-constrained parameter space 
the iE -model will leave slightly different signatures on the 
perturbation evolutions from those left by the standard 
ACDM paradigm in general relativity, and so cosmolog- 
ical data on cosmic microwave background (CMB) and 
matter power spectra might place some constraint on the 
parameter space. 



A. The Constrained Parameter Space 

In this subsection we briefly summarize the constraints 
on the M -model described in Eq. It has been well 

known that in the weak-field, slow-motion, limit and in 
the background cosmology the iE -model displays track- 
ing behavior. For the former environment, it can be 
shown that in the presence of the iE -field the observed 
effective gravitational constant Go is a rescaling of the 
bare one Gn [IOj dl by (note that here our as have 
different signs from those in p^ ) 



Correspondingly, we define — SttGo and Koo — SttGoo 
to be used below. Note that the rescaled Goo is gener- 
ally not equal to Gn and as such the background cos- 
mic expansion rate will be different from that in stan- 
dard GR. However, we note that the numerical value 
of gravitational constant we find in the textbooks and 
use in the numerical calculations is not Gn but rather 
the locally- measured Gq. It is possible to obtain lim- 
its on non-local values of G by considering the primor- 
dial nucleosynthesis of light elements (see for example 
Ref. [IH). If Ci4 = —a then we have Goo = Gq, which in- 
dicates that the background cosmological dynamics will 
be exactly the same as assuming standard GR [l^, [20[ ; 
otherwise we can use constraints from primordial nucle- 
osynthesis on the value of gravitational constant to show 
that IGoo/Go - 1| < O(O.l) [13]. However, note that 
the particle-horizon size at the epoch of neutron-proton 
freeze-out (i ~ 1 s) , which is most sensitive to variations 
in the value of G, is only ~ 10^° cm and this causally 
linked region expands in size by about a factor of 10^*^ 
by the present-day to a size ^ 10^° cm ~ 32 pc which is 
a sub-galactic scale but there has then been local grav- 
itational collapse by a factor of 10^. Such collapse may 
however also affect the local value of the gravitational 
constant (2^ , but we will not consider this in the present 
work. 

There are also constraints from the observations of pa- 
rameteriz ed p ost-Newtonian (PPN) parameters [1^ . It is 
shown in [l^] that for all the PPN parameters to coincide 
with those in GR (otherwise the parameters may need to 
be fine-tuned) one reduces the full four-parameter space 
of the model to a two-parameter subspace characterized 
by 



C2 



C4 



C1C3 



3ci 



Cl 



(37) 



Gn 



5C14 



Gn, 



(35) 



In addition, iE -theories contain five gravitational and iE 
-wave modes, which include the two usual spin-2 gravi- 
tational waves, and three additional modes: two spin-1 
transverse iE -gravity waves and one spin-0 longitudi- 
nal iE -gravity wave. The speeds of the three additional 
modes are generally not equal to 1. It has been shown in 
Ref. [l3] that if these speeds are less than 1 then the high- 
energy particles will produce Cerenkov radiation when 
passing through vacuum, which imposes stringent con- 
straints on the model. However, as suggested in p^ . 
these constraints do not apply if these modes propagate 
superluminally. The requirement that the additional iE - 
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gravity waves do not propagate subluminally further lim- 
its the parameter space to 



1< Ci3 < 0, 



Cl3 



3(1 + Ci3) 



< ci - C3 < 0. (38) 



Finally, when the above constraints Eq. ([58)1 are satis- 
fied, the positive energy requirement and the stability of 
additional wave modes fld\ also hold. In addition, we will 
have ci4 = — a so that the Big Bang nucleosynthesis con- 
straint does not apply. Thus, we can see that, even after 
using all the current constraints, there is still a large pa- 
rameter space remaining for the model. In what follows 
we shall ask whether linear-order cosmological observa- 
tions such as the CMB and the form of the matter power 
spectrum could reduce this parameter space further, and 
as we will show, the answer is positive. However, the 
modifications are small and depend weakly on the model 
parameters, which mean that the data on CMB and mat- 
ter power spectra cannot give very stringent limits on the 
parameter space. 



B. Linearly Perturbed Equations 

In this subsection we consider the perturbation evolu- 
tions of our JE -model. For generality, we will derive the 
equations for arbitrary choices of parameters and only 
later confine ourselves to the parameter space described 
in Eq. (pSj) . Besides, since the presence of the JE -field 
in general will modify the cosmology at all times, we will 
also investigate its effects during the inflationary era as 
in [l^; after that we will turn to its effects on late time 
cosmology. 

Following [l^, we shall make the following harmonic 
expansions of our perturbation variables 

^-^ a ^-^ a 

k k 
k k 

^ — ' ^ — ' a 

k k 

^-^ ^-^ a 

k k 
sea = X! ^Qa, Va^Yl —'lQa 



E 



(39) 



in which is the eigenfunction of the comoving spatial 
Laplacian a^V^ satisfying 



and QlQl, are given by Ql = aQ"" , Ql, = {aQ%- 
Note that ae is dimensionless. 

In terms of the above harmonic expansion coefficients, 
Eqs. dm m m [H mi US) can be rewritten as [3 



-fc2(cr- Z) 



nqa 
1 



Kfl^ [/c(n + x) + mq] 



(40) 
(41) 
(42) 



k{(T' + n(T) = fc2((/)-h A) - -Kfl^n, 

k^{(j)' +n(f>) = ^Ka^[k{p + p)a + kq-n' ~nn],{A3) 



k^T] = nXa^ - 2knZ, 
kri' = -Kqa^ - 2knA 



(44) 
(45) 



where Ti — a' /a — ^ad and a prime denotes the deriva- 
tive with respect to the conformal time r {adr = dt). 
Also, Eq. ([M)) and the spatial derivative of Eq. (^5)1 be- 
come 

q' + mq+{p + p)kA~kXP + \kJl ^ 0, (46) 
X' + 'ih'{p + p)+m{X + XP) + kq = (47) 

Recall that we shall always neglect the superscript *°* 
for the total dynamical quantities and add appropriate 
superscripts for individual matter species. Note that 



h' 



\kZ ~ HA 



(48) 



and that the k appearing above is the bare (not neces- 
sarily the measured) one. Furthermore, for convenience 
we define the frame-independent (FI) variables [l^ 



q+ {p + p)a, 
ee + a, 



(49) 
(50) 



^ is FI according to Eq. ([55)1 because we know that the 
anisotropic pressure tensor TTat is frame invariant. Hence, 
it follows from Eq. (fTS]) , that q is also FI up to first order 
in perturbation. In the zero-shear frame (the Newtonian 
gauge), we have simply q = q and ae = ae. 

Before presenting the evolution equations, we first 
write the dynamical quantities of the JE -field in terms of 
the harmonic coefficients introduced above: 

nX^a^ = akn{Z + a3) 



KXP^'^a^ 



nq^a^ 



-ci4 [fc^A + fc(ae' + Hffi)] , (51) 
-^fc[(Z + a3)' + 2H(Z + ae)] 



(52) 
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-ak'^{Z + £e) - -cisk'^ as, (53) 



nH^a^ = ci3fc(ae' + 2Hffi). 



(54) 



In these expressions we have used both ae and ae because 
not all these quantities are FI. However it can be shown 
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that the two quantities 
P' 



1 



P 



X 



p' J p' 
ci4 (fcae' + kHs; - 

77' 



3 3 ori^ 



(55) 



where the superscript means the fluid matter. We see 
that on large scales these attributes of the JE -field track 
those for other matter species in the universe. As the 
combination X^'^ — yX-^ determines the type of pertur- 
bations (for example, the perturbation is adiabatic if the 
combination is equal to zero), this indicates that the M - 
field will not alter the type of the scalar perturbation pro- 
duced in the inflationary era. Note that the above track- 
ing behaviors are the same as that in the background 
cosmology, i.e., rescaling the gravitational constant by a 
same factor. 



and 
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Zi o o 



(56) 



Now we can proceed to derive the evolution equations 
for our M -model. The propagation equation for the M 
-field Eq. ([^5]) . in terms of the FI variables, becomes 



which will be used in the derivations, are FI, as they are 
expressed in terms of FI variables only. Note that in the 

2 

above we have defined $ = — ffj-II for convenience, 
where $ is the Newtonian gravitational potential, and 



^ = ]^na^{p + p) 



1 



+- 



(ci4 - ci3a)e" - a(l + ci3)(W' - H'^)e (59) 



a + 2c 



13,2 



k^e- {a + cu)k{a<i>y - -a{Knfa^)' = 



We now investigate in detail Eqs. ((55l [56]) . On large 
scales, where kT <^ 1, the terms involving k (these in- 
clude the n term in Eq. ([55)1 ) can be safely disregarded, 
and as a result we have 



nq 



1 i^H 



(58) 



in which we have changed the variable to e = aa; for sim- 
plicity. Taking the time derivative of Eq. (|43p , adding to 
it (4-f 3^)7i times the same equation, and using Eq. (|46|) 
to eliminate the q' term, we obtain the following equation 
second-order differential equation for 



1 2 

-Ka 
2 



2H - — U> 



n" 



5 - — nil' + 2n'n + 6 - 2— n^n 



HE 



;„2 



3HS 



(60) 



Eqs. ([59l [60]) are the evolution equations for the coupled 
e — $ system that we are looking for (note that 11 — 
YiJE -q/ where 11-'' is the fluid matter anisotropic stress 
and XI"^ can also be expressed in terms of e and its time 
derivatives by virtue of Eq. ([51)1 ). They are not closed 

if XP'f — ^X^ and II-^ are unknown and to know these 
quantitie we would need to know the matter content of 
the universe. 



C. The Primordial Power Spectra 



In the analysis above, we have mentioned that the pres- 
ence of the jE -field does not affect the form of the scalar 
perturbation produced during inflation. But we also need 
to know whether other features of the inflationary power 
spectrum, such as the spectral index and the amplitude, 
are modified by the IR -field, as compared with the pre- 
dictions in standard GR. Here, we will investigated this 
issue (see [l^ for a study in the absence of the C4 term in 
the jE -Lagrangian) by considering a single-field model 
of inflation 24] in the presence of the -field. 
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During the inflationary epoch a scalar inflaton field ip 
slowly rolls along its potential and has an almost con- 
stant energy density which drives an almost exponen- 
tial expansion of the universal scale factor. The comov- 
ing Hubble length (the horizon) decreases with time so 
that in this process the quantum vacuum fluctuations 
of the inflaton fleld ip on the scales of interest to us 
leave the horizon (their scales become larger than the 
horizon). The curvature perturbations they generate re- 
main constant during their subsequent super-horizon evo- 
lution, until these scales eventually reenter the horizon 
long after the inflationary period has ended. During the 
radiation-dominated era when these modes stay outside 
the horizon, the metric perturbation $ becomes a con- 
stant, which drives the density perturbations of different 
matter species, and leads to the observed CMB and mat- 
ter power spectra after horizon re-entry. 

There are no couplings between the scalar field and 
gravitational or iE -fields, so we can write down its dy- 
namical quantities for the scalar field Lp as 

- \p^-V{p,), 

- 0. (61) 
Making the following harmonic expansion for Va<y5 

a 

k 

I 



it is easy to get 

= ^{p'x' + ^"A + a'V^x), (63) 

XP^^ - ^{p'x' + ^'^A-a'V^x), (64) 

q'f = -^kif'x (65) 

where = dV{ip)/dip. Then parallel to Eqs. ([HSl [HI) 
we have, following the standard procedure. 




- ^ (l + W'x' - ¥'"1' + a'V^x + 3Wx]m 

and 

= nkp'x (67) 
where we have defined the FI variable 

X = X + j'J- (68) 

Substituting Eqs. jMl ET]) into Eq. ([601), and using 
Eq. (|4ip to eliminate the term proportional to p'x' — 
p'(S? + a^VipX + STii/s'x? we arrive at the following equa- 
tion 



1 



1 + -CM I fc2$ 



1 



1 1/ 

-Ci4 - - {C1 + C2 + C3) 



Cl3 



k , 

— e 

a 



1 , , S' fc 1 / e- 

(Ci +C2+ C3) — -e - -Ci3 ( — 



S- . (69) 
a J 



Eqs. ([591 [69|) that form a closed set of evolution equations 
for the coupled iE -inflaton system; they are a general- 
ization of the results presented in jl6| , and from them we 
can perform our analysis of the observational effects. 

Let us look first at the large-scale evolution of 4>. In 
this limit fcr ^ 1 so the last term on the left-hand side, 
and the first two terms on the right hand side, of Eq. ([69]) 
can be dropped. Notice that during the inflationary era 
rC^ -n' ^ \Ka'^{p + p) ~ and Eq. ([59]) become 



(ci4 - ci3a)e 



(a -I- ci4)fc(a$)'. 



(70) 



Substituting this into Eq. ([69[) one could see that 



k 



Cl3 



a 



a 



^ a 



n $ 



so that we finally have 

+ (21-1 - ^ ) $' + ( 2n' - ) $ = (71) 



on super-horizon scales. Note that one should not simply 
set e" = in this limit as was done in Ref . [l^ , although 
that leads to the same result [28|. This equation is still 
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valid outside the inflationary epoch within the parameter 
space given by Eq. ([55]) . 

The solution to Eq. ([7T|) is given by 



$ = £> 



1-^1 a'dr 



(72) 



where D is a constant. Long after leaving the horizon, 
the potential $ and the scalar perturbation x have time 
to evolve, and after inflation x ceases to exist. How- 
ever, as mentioned above, during the radiation era on 
super-horizon scales $ finally becomes a constant which 
is related to D by 



$ = -D. 



(73) 



Once D is known we could set the initial conditions of <& 
in the radiation era for the subsequent evolutions. The 
quantity D could be fixed by matching to the inflaton 
perturbation x at the time of horizon exit (fc = aH — Ti.) 
as follows: substituting Eqs. ([551 [57]) into Eq. we get 



fc(a$)' 



1 



1 - ia 
-C13 W- 



^kaKifi'x - ^(ci +C2 + C3)fc^e 
f {H' - n^)e] . (74) 



Now we begin to confine ourselves within the parame- 
ter space Eq. ([57]) . In that case Ci4 = — a = (ci +C3)(ci — 
C3)/ci and the JE -field EOM Eq. ^ becomes 

a{l + cu)[e" + in' -n^)e] = {ci + C2 + C3)kh.i75) 

So, from these two equations we find that 

1 1 



(a$)' = 



1 



— -anip'x. 



On the other hand, from Eq. (|72p one can write 

1 



21 



(76) 



(77) 



Obviously, matching these two expressions gives the value 
of £> 



D 



—x 

ip' 



(78) 



where the three variables Ti., tp' and x are all evaluated at 
the horizon exit. As a result, the initial power spectrum 
for = A;3($2)/27r2 where (• • • ) means the ensemble 

average, is given by 



4 7/2 



4W^ 

9 (^'2 



(79) 



in which we have used the relation = (x^) /2tt'^ = 
{H/2itY HI] (as discussed in Ref. this relation is 



not affected by the jE -field up to first order because the 
M -field is not coupled to the inflaton) and the Hubble 
expansion rate H — a /a is also evaluated at the time of 
horizon exit. 

Note that the $ power spectrum Eq. ([79)) has exactly 
the same form as in standard GR. To compare their mag- 
nitudes, let us write 



{Hyp''- 



{Hyp'- 



in which we have used the fact that at the time of horizon 
exit Ti — k should be the same in the two models. If 
we further assume the same inflation potential, then [20| 
because the background expansion in our jEther model 
is indistinguishable from that of GR, we have 



(80) 



that is, the power spectra in these two models are exactly 
the same, in both the shape and the magnitude. It should 
be noticed that if the background evolution, of the 
jE-model during the inflation is different from GR, then 
the evolutions of p will also be different because of the 
scalar field equation of motion; as a result the slow-roll 
parameters at the horizon crossing are generally different 
in these two models and so will be the spectral indices. 

What is the fate of the vector and tensor modes of the 
jE -field perturbation? As shown in [l^ , the vector mode 
is also sourceless unless exotic matter such as cosmologi- 
cal defects exist. Hence, it should decay similarly all the 
way up to the present and leave no traces in the CMB 
observables. As for the tensor modes, from Eqs. ([29l - [33)) 
one can see that the contribution of the M -field to the 
tensor modes lies only in kH^ and is independent of the 
model parameter C4. Consequently, the power spectrum 
of the tensor perturbations given in [l^ 



T'h = k(1 + ci3 



,1/2 



H^ 
27r2 



(81) 



is still valid here (note that what appears in this formula 
is the true bare gravitational constant). Following the 
argument above again, we can show that 



-pGR 
' h 



I+2C14 



(1 + C13)1/^. 



(82) 



According to Eq. ([38)1 . C13 < and C14 < 0, so the ampli- 
tude of the gravitational-wave spectrum predicted in our 
jE -model is smaller than that arising in standard GR. 
Although this may provide a discrimination between the 
two models, the difhculty in observing the tensor spec- 
trum will be a great hurdle for the use of this discrimi- 
nator: if our Universe turns out to be described by the 
jE -model, then this spectrum is even more difficult to 
observe than in GR. 
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whose solution is 




FIG. 1: (Color online) The evolution of the M -field pertur- 
bation ae versus the cosmic scale factor a. We have shown it 
for 3 difi'erent scales k = 0.1,0.01,0.001 Mpc~^ respectively, 
as indicated beside the curves. The model parameters are 
a = 0.2, ci3 = —0.2, which satisfy the constraints Eq. (|38}. 



D. The M -Effects on Late-time Cosmology 

We have seen in the above that in our M -model the pri- 
mordial scalar-mode power spectrum is exactly the same 
as assuming standard GR and using the measured value 
in calculations. The next question is if possible devia- 
tions from GR will emerge in the subsequent evolutions of 
the perturbations. Interestingly, we find that the answer 
is yes. The reason is that, although the iE -perturbation 
is sourceless during the inflationary era, in subsequent 
stages where the fluid anisotropic stresses are nonzero it 
becomes sourced. To see the consequences explicitly, we 
now write the jE -field EOM within the parameter space 
Eq. dSH) as 



a [(1 + ci3)fc(ae' + + kH^o^ 



a + 2ci3 



3 

During the inflationary epoch and H'- 

thus the solution to the above equation is 

3s{t) ~ — ^- sin (csfcr) 
a 

where we have defined ci 



ae. 



~ 0, 



ffigcn(T) = ^ [cos(csfcr)csA;T - sin(csfcr)] 
B2 

H — ^ [sm{cskT)cskT + cos{cskT)] (85) 



where Bi^2 are integration constants and a subscript gen 
means the general solution. We can see that in the limit 
Cgkr < 1 afegcn decays as ajgcn ~ l/r^ and when Cskr ^ 1 
it decreases as cKgon ~ -Bi cos(csfcr)/T + B2 smicskr) / t . 
Meanwhile there is a particular solution of ac which in- 
volves a weighted integration of [nlif a?)' over time. Deep 
into the radiation-dominated epoch, the decaying general 
solutions of ae have become negligible and the growing 
particular solution is still tiny, so we can reasonably take 
the initial conditions of = in our nu- 

merical calculations. We have checked that, for the scales 
we are interested in, this choice of initial conditions is ro- 
bust and the numerical results are not sensitive to (not 
dramatically) different initial conditions, and in Fig. [1] we 
have depicted the time evolution of is at different scales 
(or different k) for the model a = — C13 — 0.2. It can be 
seen there that, at early times jb remains close to zero; 
later it grows as {kH^ a?)' deviates significantly from 0, 
and finally when {kH^ a?)' becomes tiny in the matter- 
dominated era it undergoes a (oscillatory) decay again. 



Now what kind of signatures will this behavior of im- 
print on the perturbation evolution? In order to answer 
this question, we find it useful to rewrite the total den- 
sity perturbation, pressure perturbation, heat flux and 
anisotropic stress (with the contributions from the M- 
fleld included) as follows: 



Klltra^ -f KllntrO^ 



(86) 
(87) 
(88) 
(89) 



a 

Cl3 



COS {cgkr) (83) where the quantities with subscript tr 



1TTY-, u T which is posi- 

3(1 + Ci3)(ci-C3) 1^ 

tive according to Eq. ((38|) : ^1,2 are integration constants. 
Meanwhile, during this epoch the scale factor a under- 
goes exponential growth, which means that ae decays ex- 
ponentially and its initial value is washed out soon. 

In the radiation dominated epoch, a (x t and n-'^ = 
W +n'' ^ (here subscripts J, u denotes photon and 
neutrino respectively) . Now the homogeneous part of the 
M -field EOM becomes 



2 



KQtra 
Klltra^ — 



K^XP-W^, 



(90) 
(91) 
(92) 
(93) 







(84) 



are the tracking parts because they are exact rescalings 
of the corresponding quantities for the fiuid matter [we 
remind that = k/(1 — ^a)], while those with subscript 
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are the non-tracking parts expressed as 



Kqntra 



a 



+ r—{l + ci3)k{S;' + 2n&), (94) 

1-2" 

-^-^(l + ci3)fc(ae' + 2Wae), (95) 
o 1 — T^a 



la + 2ci3 2~ 
3 1 - ia 



(96) 



ci3fc(ae' + 2'HSi) — -q:k^ n-'^a^; (97) 



only the measured value appears in the final expres- 
sions. It can be easily checked that the tracking and 
non-tracking parts satisfy the energy momentum conser- 
vation equation separately by utilizing the M -field EOM. 
In particular, the non-tracking part effectively represents 
a purely perturbed contribution to the total energy mo- 
mentum tensor that has no background counterparts. We 
also note that kX^^^^o? = —^KX^t^a? . 

Taking the time derivative of Eq. ([44]) , using Eq. ([45]) 
to eliminate the r/ term and cancelling the X' appearing 
there with Eq. (|47p. one arrives at the following evolution 
equation of Z: 



k{Z' +nZ) + ^K(X + 3XP)a^ = 



or equally 



k{Z' + nZ) + ^K,4X^ +3XP'f)a^ = 0, 



(98) 



(99) 



that is, the non-tracking part does not contribute to 
the evolution of 2^. In this subsection we work in the 
CDM frame where the 4- acceleration A — 0, and thus 
^CDM = —kZ [18] (where A = X/p is the density 
contrast) which means the growth of matter perturba- 
tion is not affected by the non-tracking part. However, 
since nXntia^, which is nonzero at early times, appears in 
Eq. ([33]), the values of K^,X^a^ and Z in this equation are 
slightly modified that change the subsequent evolution of 
Z a little. These can be seen in Fig. [31 where we have 
displayed the matter power spectrum of the M -model 
for different choices of parameters. Obviously the spec- 
trum depends weakly on the parameters; furthermore, its 
shape is essentially the same as that in ACDM model so 
that the parameters will degenerate with the bias relat- 
ing the matter power spectrum to the observed spectrum 
of galaxies. In addition to that, there is also the degener- 
acy with respect to the neutrino masses, since the latter 
affects its anisotropic stress and the time at which they 
become non-relativistic. Taken together this indicates 
that the data for P{k) will not place very stringent con- 
straints on the M -model. 

From Eqs. ([3I]-[3S|) one can see that are also 

influenced by the non-tracking part, and because these 



t 




k (Mpc ') 



FIG. 2: (Color online) The matter power spectrum of the 
M -model considered in this work. From top to bottom the 
curves correspond to (0,013) = (0.5, —0.5), (0.2, —0.2) and 
(0, 0) (which is just ACDM) respectively. The other parame- 
ters common for all curves are Qth^ = 0.0223, flch^ = 0.1054 
and h = 0.732 where Qb,c are the fractional energy density 
of baryons and CDM, and h = Ho/{WO km/s/Mpc). The 
differences among these spectra are small although their pa- 
rameters are quite different. 



variable determine the CMB power spectrum, we could 
expect the latter to be modified by the presence of the M 
-field as well. This is correct, as shown Fig. [31 where we 
have plotted the CMB spectra for different choices of the 
model. Obviously the CMB power spectrum is changed 
by the M -field, though this effect is not very strong and 
so it might be difficult to place stringent limits on the 
parameter space of the theory. 

There are two points that need to be noticed. Firstly, 
in this work we are assuming that the recent cosmic 
acceleration is caused by a cosmological constant, the 
anisotropic stress of which is zero. If in contrast the dark 
energy has a nonzero anisotropic stress (as in [25l[26ll27|l 
then the M -field perturbation ae probably will not decay 
at late times, and this might lead to further modifica- 
tions of the CMB power spectrum. Secondly, when f{K) 
is not simply a linear function of K, the JE -field will be 
sourced by the evolution of the gravitational potential as 
well, and its effects will be more intricate and interest- 
ing. Such new features may appear in more complicated 
models like TeVeS 0, Q and JE -fields with non-canonical 
kinetic terms [6^. 

In total, these results indicate that within the param- 
eter space Eq. ([571) the background cosmology of our M 
-model is the same as in GR, while the CMB and matter 
power spectra differ slightly from the predications of GR. 
In principle, these features could be used to constrain the 
parameters a and C13 of the present model. This will 
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FIG. 3: (Color online) The CMB power spectrum of the M 
-model considered in this work. At i ~ 10 from top to bottom 
the curves correspond to (a, C13) = (0.5, —0.5), (0.2, —0.2) and 
(0, 0) (which is just ACDM) respectively. The other parame- 
ters common for all curves are Qbh^ = 0.0223, Qch^ = 0.1054 
and h — 0.732 where Qb,c are the fractional energy density of 
baryons and CDM, and h = //o/(100 km/s/Mpc). 

generally involve a full search of the parameter space us- 
ing for example the Markov Chain Monte Carlo method, 
which is beyond the scope of this work. However, as we 
discussed above, the constraints from the linear spectra 
may not be very stringent anyway. 



IV. DISCUSSION AND CONCLUSION 

In this paper we have studied the cosmology of the 
Einstein- iE ther theory. After presenting the general field 
equations for such theories in the CGI formalism, we fo- 
cussed on a specific class of models described in Eq. (l34l) , 
and confined ourselves to the parameter space of models 
described by Eq. ([55)1 which pass the PPN and Cerenkov 
constraints. This parameter space is known to have the 
same locally- and cosmologically-felt gravitational con- 
stants (kq = Koo, which are different from the true bare 
k) and this tracking behavior indicates that we can con- 
sider its background expansion just by using the mea- 
sured value Ki, = kq and ignoring the presence of the M 
-field (since the only effect of the iEther field is to track 
other matter species, the model in Eq. clearly can- 
not explain dark energy. More general models, such as 
the f{K) proposed in @ and our Sec.|lTl may serve this 
purpose). 

We find that it is a general feature that the tracking be- 
havior not only occurs at the background level but also 
at the linear order in perturbation theory. For exam- 
ple, the quantity XP'^ — I^X-^ tracks that of the other 



matter species on super-horizon scales. This indicates 
that, whatever type of perturbation is generated during 
inflation, the presence of M -field will not alter it. In 
particular, in the single-field inflation model we consider, 
no isocurvature perturbation is produced. This is an im- 
portant characteristic, and it would be interesting to see 
whether similar behavior occurs for general, higher-order 
choices of f{K). 

We generalized the analysis of primordial power spec- 
tra for the JE -theory to our model. For the param- 
eter ranges which satisfy local gravity bounds, we find 
that the evolution of the large-scale gravitational poten- 
tial, <&, is unmodified as compared with that in GR, and 
show that the primordial power spectrum of also has 
the same form as in GR [c.f. Eq. ([75)1 ]. If we assume that 
the bare gravitational constants in GR and in the JE - 
model are the same, then the magnitudes of the spectra 
are different in these two models. However, contrary to 
the discussion in [l6l | , we argue that we do not know the 
true bare k, but only know the measured value, k^, — kq. 
In both GR and our JE -model, kq is equal to the cos- 
mological value Koo, while in the latter k ^ Kqo- As a 
result, we show that with the same inflationary potential 
the primordial power spectra of $ in these two models 
should have the same shape and the same magnitude. 
Meanwhile, we also find that the power spectrum of ten- 
sor perturbations in our model is smaller in magnitude 
than is predicted in GR, and so currently we cannot use 
it to distinguish between GR and the JE -model. 

For the late-time evolution of the perturbations, it is 
shown that, the IE -field perturbation is driven by the 
evolution of the anisotropic photon and neutrino stresses 
- the [K*(n^ + Il^)a^]' term in its propagation equa- 
tion - although it is sourceless during inflation. There- 
fore, it decays away exponentially during the inflation- 
ary era, grows again in the radiation-dominated epoch 
when [k^(H'' + Il'^)a^]' is significant, and finally dimin- 
ishes again, oscillating when [/{^(H'^ + Il'^)a'^y eventually 
becomes negligible. As a result, the CMB and matter 
power spectra are modified by the existence of the JE 
-field. We also remark that, depending on the nature 
of the dark energy, it is possible that the JE -field per- 
turbation will have non-trivial dynamics driven by the 
{Ki^lP^a^y term, where H^^ is the possible dark-energy 
anisotropic stress [25| . 

We also note that there recently appeared a later paper 
[2^ which investigated in details the structure formation 
of the model proposed in In their Sec. V the authors 
considered a simple power-law model f{K) = ^{—K)^ in 
which n = 1 corresponds to our Lagrangian Eq. p4p (ex- 
cept for the C4 term). They obtained criteria for there to 
be a growing mode in the evolution of the (scalar-mode) 
jE-field perturbation. However their criteria cannot be 
applied directly to the model we consider here. To see 
why, note that in Eq. (32) of [2^ the quantities \E', $ also 
contain V (their V is equivalently our ae). However, by 
their argument the V terms in ^f, $ are suppressed by 
a small quantity Fk (which in our notation is just F) 
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and are ultimately neglected from the jE-field equation 
of motion, their Eq. (43). In our model, there is no ar- 
gument that F 1 (its magnitude actually is 1) and so 
we can no longer neglect the ae terms appearing in ^P. 
When these are taken into account we obtain a different 
equation of motion, i.e., the biS in Eq. (43) of [2§| are 
different in our work, and consequently the criteria for 
the existence of growing modes are different as well. In 
addition, are the facts that we have a C4 term in our La- 
grangian and use a different parameter space, Eqs. p71 
, which also contribute to the differences between our 
results and those in [2^. As an aside, we stress that 
the key relation in our work, C14 = —a, is a consequence 
of Eq. (l37|) and has nothing to do with whether or not 
the iE-gravity waves propagate superluminally. In fact, 
we could do our calculation dropping the constraint that 
these waves propagate superluminally, and in this case 
the perturbation dynamics might place some constraints 
on the parameter space; for example, in some portion of 
the parameter space the growth of ae may become un- 
stable. In our work we choose the parameters leading to 
superluminal gravity waves simply to avoid constraints 
from Cerenkov radiation. There are still some debates 
about this kind of choice, see however [l^l for a different 
and conservative point of view. 

There are also other general differences between the 
EA model we consider here and the GEA model of ^29i] . 
First, for the static and weak field limit, the C4 term and 
the choice of parameter space Eq. (l37|l arc crucial for 
the EA model to evade PPN tests ^10]; in '2^ there is 
no C4 term, but the nonlinearity in the iEther field La- 
grangian guarantees that the local gravitational tests are 
not a problem because at high densities the modification 
is negligible. Second, as we have seen above, our con- 
clusion that the primordial power spectrum of density 
perturbation is unmodified compared with GR relies on 
the fact ci4 = — a so that the background expansion dur- 
ing infiation in our model is the same as in GR (again the 
C4 term is crucial here). If the C4 term is not included, the 
local and cosmological gravitational constants are differ- 
ent, meaning that the background evolution during in- 



flation is different from GR; as a result the background 
quantities at horizon crossing, which determine the ob- 
servables in inflationary models through the slow-roll pa- 
rameters, may be different from those in GR - this means 
that the EA model will predict different shape of primor- 
dial density spectrum other than GR 30]. For the GEA 
model considered in 29], again the nonlinear f{K) makes 
the modifications to GR during early times small enough 
to be neglected, and as such the primordial power spec- 
trum is also the same as in GR. Note however that in 
the latter case [2^ the late time evolution of vector field 
perturbation is also modified in a complex way, making 
the corresponding cosmological behaviors different from 
ours. 

So, in conclusion, although the background expansion 
of our JE -model is exactly the same as that in GR, the 
cosmological data on CMB and matter power spectra can 
be used to constrain the JE -model. However, the con- 
straint is not expected to be very stringent because the 
linear perturbation spectra depend weakly on the model 
parameters. Other considerations of the behaviour in 
strong gravitational fields, such as those studies of the 
compact stars or black holes [H, [3| would enable better 
constraints on the present model to be obtained. Alao, as 
general interests for cosmology, studies in more compli- 
cated iE-field Lagrangian as those performed in [2^, [3l| 
need to be explored thoroughly. 
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